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ANALYSIS  OF  HEAT  TRANSFER  AND  PRESSURE  DROP  FOR  A  GAS 
FLOWING  THROUGH  A  SET  OF  MULTIPIE  PARALLEL 
FIAT  PLATES  AT  HIGH  TEMIERATURES 
By  Thomas  H.  Einstein 


SUMMARY 

Equations  were  derived  representing  heat  transfer  and  pressure  drop 
for  a  gas  flowing  in  the  passages  of  a  heater  composed  of  a  series  of 
parallel  flat  plates.  The  plates  generated  heat  which  was  transferred 
to  the  flowing  gas  hy  convection.  The  relatively  high  temperature  level 
of  this  system  necessitated  the  consideration  of  heat  transfer  between 
the  plates  by  radiation. 

The  equations  were  solved  on  an  IBM  704  computer,  and  results  were 
obtained  for  hydrogen  as  the  working  fluid  for  a  series  of  cases  with  a 
gas  inlet  temperature  of  200°  R,  an  exit  temperature  of  5000°  R,  and  exit 
Mach  numbers  ranging  from  0.2  to  0.8.  The  length  of  the  heater  composed 
of  the  plates  ranged  from  2  to  4  feet,  and  the  spacing  between  the  plates 
was  varied  from  0.003  to  0.01  foot.  Most  of  the  results  were  for  a  five- 
plate  heater,  but  results  are  also  given  for  nine  plates  to  show  the 
effect  of  increasing  the  number  of  plates.  The  heat  generation  was  as¬ 
sumed  to  be  identical  for  each  plate  but  was  varied  along  the  length  of 
the  plates.  The  axial  variation  of  power  used  to  obtain  the  results  pre¬ 
sented  is  the  so-called  "2/3-coslne  variation."  The  boundaries  surround¬ 
ing  the  set  of  plates,  and  parallel  to  it,  were  assumed  adiabatic,  so 
that  all  the  power  generated  in  the  plates  went  into  heating  the  gas. 

The  results  are  presented  in  plots  of  maximum  plate  and  maximum  adi¬ 
abatic  wall  temperatures  as  functions  of  parameters  proportional  to 
f(L/D),  for  the  case  of  both  laminar  and  turbulent  flow.  Here  f  is  the 
Fanning  friction  factor  and  (L/D)  is  the  length  to  equivalent  diameter 
ratio  of  the  passages  in  the  heater.  The  pressure  drop  through  the  heat¬ 
er  is  presented  as  a  function  of  these  same  parameters,  the  exit  Mach 
number,  and  the  pressure  at  the  exit  of  the  heater. 

One  of  the  more  interesting  effects  associated  with  this  heater  ge¬ 
ometry  was  the  nonuniform  distribution  of  weight  flow  between  the  pas¬ 
sages,  despite  the  fact  that  the  heat  generation  for  adl  plates  was  the 
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same  and  that  the  spacings  between  the  plates  were  identical.  It  was 
found  that  the  magnitude  of  this  nonuniformity  of  flow  distribution  in¬ 
creased  with  decreasing  f(L/D).  Another  unusual  effect  was  that,  when 
the  f(L/D)  ratio  of  the  heater  was  decreased,  although  the  maximum  plate 
temperature  became  higher  as  expected,  the  maximum  adiabatic  wall  tem¬ 
perature  actually  decreased  substantially. 


INTRODUCTION 

The  analysis  of  heat  transfer  by  forced  convection  to  a  fluid  flow¬ 
ing  in  a  single  passage  between  a  pair  of  parallel  flat  plates  is  well 
known  for  the  cases  of  both  turbulent  and  laminar  flow,  and  is  thorough¬ 
ly  treated  in  the  literature.  However,  little  work  has  been  done  on 
treating  the  problem  of  heat  transfer  to  flow  in  a  series  of  parallel 
passages  formed  by  a  set  of  multiple  parallel  flat  plates.  The  basic  re¬ 
lations  governing  heat  transfer  from  the  plates  to  the  fluid  in  the  pas¬ 
sages  remain  unchanged  from  those  in  a  single  passage;  however,  the  prob¬ 
lem  is  now  complicated  because  there  will  be  interaction  between  the  heat 
transfer  to  one  passage  and  that  to  adjacent  passages,  resulting  in  non- 
\miform  plate  temperatures  and  nonuniform  flow  distribution  between  the 
passages.  In  addition  to  the  foregoing  complications  that  arise  in  the 
multipassage  case,  if  the  heater  formed  by  the  set  of  parallel  plates  is 
operated  at  a  high  temperature,  then  radiation  between  the  plates  becomes 
significant.  This  effect  must  also  be  considered.  It  is  therefore  the 
purpose  of  this  report  to  present  an  analysis  of  this  heat-transfer  ge¬ 
ometry,  and  tc  obtain  some  quantitative  results  at  high  temperatures, 
which  reflect  the  effect  of  interaction  between  the  passages  and  teike 
into  consideration  the  effect  of  radiation  between  the  plates. 

Before  proceeding,  one  might  justifiably  raise  the  question  of  where 
this  heat-transfer  geometry  would  find  application.  In  the  general  case, 
the  plates  would  be  heated  either  electrically,  or  by  nuclear  fission, 
the  latter  alternative  implying  the  use  of  this  geometry  as  a  fuel  ele¬ 
ment  in  a  nuclear  reactor.  Many  of  the  stationary  nuclear  reactor  power- 
plants  in  use  today  use  this  scheme  for  a  fuel  element,  and  the  applica¬ 
tion  of  this  geometry  for  a  fuel  element  is  of  interest  for  nuclear  rock¬ 
et  reactors. 

This  report,  then,  describes  an  analysis  of  heat  transfer  and  pres¬ 
sure  drop  in  a  high-temperature  heater  consisting  of  parallel  flat  plates. 
The  analysis  considers  the  effects  of  radiation  between  the  plates  and 
allows  for  both  laminar  and  turbulent  flow  of  the  working  fluid.  Provi¬ 
sions  are  also  made  for  adjusting  the  division  of  fluid  flow  between  the 
passages  such  that  the  static-pressure  drops  in  al  1  the  passages  will  be 
equalized. 
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Specific  results  will  be  presented  using  hydrogen  as  the  gas  flowing 
through  the  heater.  Although  hydrogen  partially  dissociates,  resulting 
in  marked  Increases  in  heat-transfer  coefficient  at  the  higher  tempera- 
t\ires  for  which  the  results  are  presented,  and  is  thus  in  this  respect  not 
typical  of  the  behavior  of  most  gases  at  more  moderate  temperatures,  this 
gas  was  chosen  for  the  presentation  of  results  because  of  the  high  level 
of  current  Interest  in  it  as  a  propellant  for  space  vehicles. 


SYMBOIS 

Op  specific  heat  of  gas,  Btu/lb-°R 

D  equivalent  diameter  of  passsiges,  ft 

^  radiation  form  factor,  Btu/ft^-sec-°R^ 

f  Fanning  friction  factor,  dimensionless 

G  mean  mass  flow  of  gas,  lb/{sec-ft2) 

g  acceleration  of  gravity,  32.2  ft/sec^ 

h  convective  heat-transfer  coefficient,  Btu/sec-ft^-^R 

J  totSLl  n\ffliber  of  plates  in  heater 

k  thermed  conductivity  of  gas,  Btu/sec-ft-^R 
L  total  length  of  plates,  ft 

M  Mach  number,  dimensionless 

N  number  of  passages  in  one  symmetrical  half  of  heater 

Nu  Nusselt  number,  dimensionless 

P  overadl  static  pressure  drop,  Ib/ft^ 

Pr  Prandtl  number,  dimensionless 

p  loced  static  pressure,  Ib/ft^ 

q  heat  flux,  Btu/sec-ft^ 

R  gas  constant,  ft/°R 

Re  Reynolds  n\]mber,  dimensionless 
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s  plate  spacing,  ft 

T  surface  or  plate  temperature,  °R 
t  total,  gas  temperature,  °R 

ts  static  gas  temperature,  °R 

V  gas  velocity,  ft/sec 

w  passage  weight  flow,  Ih/sec 

X  distance  from  inlet,  ft 

T  ratio  of  specific  heats,  dimensionless 
€  radiation  emissivity,  dimensionless 

(i  viscosity  of  gas,  Ih/ft-sec 

p  desnity  of  gas,  Ib/ft^ 

o  Stefan-Boltzmann  constant,  4.75X10"^^  Btu/sec-ft^-°R^ 

Subscripts : 
b  bulk  conditions 

f  film  conditions 

g  generated 

i  axial  increment  number 

in  inlet  conditions 

J  trial  nimber  in  iteration  for  pressiire  drop  equalization 

Z  trial  number  in  iterative  solution  of  plate  heat  balance  equations 

n  p^-ate  or  passage  number 

o  exit  conditions 

X  axial  location 

Superscripts: 

_  mean  value 

'  from  interior  side  of  a  plate 
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ANALYSIS 

A  sketch  of  the  parallel  plate  heater  under  consideration  is  shown 
in  figure  1.  The  analysis  of  its  operation  will  he  divided  into  three 
distinct  sections.  First,  the  analysis  of  combined  forced  convection  eind 
radiation  heat  transfer  from  and  between  the  plates  will  be  presented, 
and  the  many  simplifying  assumptions  that  are  made  in  this  area  will  be 
discussed.  The  second  phase  of  the  analysis  to  be  covered  is  the  calcu¬ 
lation  of  the  combined  friction  and  momentum  pressure  drop  in  each  pas¬ 
sage.  Finally,  the  iterative  method  of  readjusting  the  assijmed  weight 
flows  in  each  passage  to  cause  eqvialization  of  the  pressure  drops  in  all 
the  passages  is  Illustrated. 


Heat  Transfer 

Because  of  the  allowance  in  the  analysis  for  variable  fluid  proper¬ 
ties,  variable  heat-transfer  coefficient,  arbitrary  variations  in  power 
distribution,  and  the  inclusion  of  radiation  in  this  heat-transfer  prob¬ 
lem,  a  closed-form  solution  is  out  of  the  question.  Therefore,  the  solu¬ 
tion  must  be  obtained  numerically  by  breaking  the  heater  up  into  finite 
segments  and  solving  the  heat-balance  equations  separately  for  each  seg¬ 
ment.  The  heater  was  therefore  broken  up  into  100  axial  segments,  and 
the  continuous  variation  of  plate  surface  temperature  and  gas  temperature 
along  the  length  of  the  heater  is  then  approximated  by  the  temperature 
distribution  of  these  100  segments.  The  segments  are  considered  small 
enough  that  the  plate  and  gas  temperatures  of  each  segment  may  be  con¬ 
sidered  to  be  constant  over  the  length  of  that  segment.  The  number  100 
was  chosen  arbitrarily  as  large  enough  to  ensure  a  sufficiently  accurate 
approximation  to  the  actual  temperature  distribution  over  the  length  of 
the  heater,  and  yet  not  so  large  as  to  require  an  excessively  long  time 
for  the  computation  of  the  solution.  Each  segment  of  length  is  then  com¬ 
posed  of  J  plates  and  J  +  1  passages.  In  order  to  simplify  the  com¬ 
putations  associated  with  the  analysis,  the  number  of  plates  J  was  ar¬ 
bitrarily  taken  to  be  odd,  and  the  arrangement  of  plate  spacings  and 
power  generations,  though  not  necesseu-ily  all  identical,  was  assumed  to 
be  symmetrical  about  the  center  plate. 

Therefore,  it  becomes  necessary  to  apply  the  analysis  to  only  one- 
half  of  the  segment,  from  the  center  plate  to  the  outer  wall.  The  ar¬ 
rangement  of  plates  in  each  segment  that  is  considered  in  the  analysis 
is  illustrated  in  figure  2.  The  heat  transfer  in  each  segment  is  by  radi¬ 
ation  between  the  plates,  and  by  convection  from  the  plates  to  the  gas 
flowing  in  the  adjacent  passages. 

For  radiation; 


^ad.  “  ■  ^n-l) 


(1) 
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The  coefficient  ^  is  a  radiation  form  factor  and  is  a  function  of  the 
plate  emissivity  and  the  geometry  under  consideration.  Its  application 
is  subject  to  several  assumptions  that  will  be  discussed  in  a  later  para¬ 
graph  . 

For  convection: 

‘Iconv  ~  ^("^n  “  '^'n)  heat  transfer  to  outer  passage 

>  (2) 

'Iconv  “  Ni^^n  "  '^i-l)  heat  transfer  to  inner  passage J 

Here  the  coefficient  h  is  the  convection  heat-transfer  coefficient.  It 
is  a  function  of  the  gas  properties,  the  spacing  between  plates,  gas  flow 
rate,  and  whether  or  not  the  flow  is  turbulent. 

Since  a  given  plate  transfers  heat  by  radiation  to  the  two  adjacent 
plates  on  either  side  of  it,  and  by  convection  to  the  gas  in  the  two  pas¬ 
sages  adjacent  to  it,  the  heat  balance  on  a  single  plate  may  be  written 
by  combining  equations  of  the  form  (l)  and  (2)  and  referring  to  figure  2. 
Such  a  heat  balance  for  the  n^^  plate  in  one  of  the  axial  Increments  is 
then: 

‘Ign  +  ^n-l(Ttl  -  -  T^+l)  -  hA-lCTn  -  t^.i) 

-  hn(Tn  -  tn)  =  0,  n  =  1,  2,  .  .  .  ,  N  +  1  (3) 

Here  N  is  the  number  of  plates  in  the  symmetrical  half  under  con¬ 
sideration  including  the  center  plate,  and  is  the  heat  generated 

in  the  n^^  plate  of  the  segment.  Since  the  centerline  of  symmetry 
splits  the  center  plate  in  half,  heat  transfer  from  only  one  side  of  the 
center  plate  is  considered  in  equation  (3),  and  the  value  of  q^^,  the 

heat  generation  term  for  the  center  plate  in  equation  (3),  is  only  one- 
half  the  total  heat  actually  generated  in  the  center  plate.  Since  heat 
transfer  from  only  one  side  of  the  center  plate  is  considered,  all  vari¬ 
ables  having  zero  for  a  subscript  (those  whose  subscripts  are  n-1,  for 
the  case  n  »  l)  should  be  considered  to  be  identically  zero.  The  sub¬ 
script  N+1  refers  to  the  outer  weLLl,  and  here  similarly  sLll  variables 
having  a  subscript  of  N+2  should  be  considered  to  vanish  identically. 

It  is  assumed  in  finding  the  solution  to  the  set  of  equations  (3) 
that  the  only  unknowns  are  the  plate  temperatures  T^j  the  other  veuri- 
ables  in  equation  (3)  are  either  known  explicitly  or  are  computable  from 
other  relations.  The  problem  of  finding  the  is  then  one  of  solving 

the  N+1  nonlinear  algebraic  equations  (eqs.  (3))  simultaneously  for 
Tl,  T2,  .  .  .  ,  Tn+1-  Tte  method  of  solution  employed  is  presented  In 
the  appendix. 

To  obtain  the  gas  temperature  rise  over  each  increment,  a  set  of 
heat  balance  equations  may  be  written  for  the  gas  flowing  in  each  passage 
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of  the  axial  segment: 

WCp  Atjj  =  -  tjj)  +  h)^^j_(Tfj^2  -  tj^)  n  =  1^  7,,  .  .  N  (4) 

Thus,  the  temperature  of  the  gas  flowing  in  each  passage  may  he  found  at 
each  increment,  once  the  plate  temperatures  are  computed  from  equation 
(3).  A  few  words  of  clarification  are  appropriate  here.  The  gas  temper¬ 
ature  in  each  passage  of  an  increment  will  of  course  continually  increase 
over  the  length  of  that  increment.  Therefore,  the  used  in  equations 

(3)  and  (4)  actually  represent  mean  values  of  t^  .over  the  Increment, 
whereas  the  At^  in  equation  (4)  represent  the  temperature  increases 
over  that  increment;  in  other  words,  At^  i  =  -  t^  i,  where  tn,  i 

is  now  the  gas  temperature  in  the  n^^  passage  entering  the  1'^^  incre¬ 
ment.  To  simplify  the  use  of  equation  (4)  in  determining  the  mean  value 
tn  in  equation  (3),  for  the  next  increment,  the  following  assumption  is 
made: 


Atn,  i  -  i+1  -  tn,  i  =  tn,  i+i  -  tn,  i  (5) 

This  is  a  good  assumption  if  the  size  of  the  increments  is  relatively 
small  (heater  broken  into  a  large  number  of  increments).  To  summarize 
the  computational  procedure,  at  each  increment  the  set  of  equations  (3) 
is  solved  simultaneously  for  the  Tnj  then  the  mean  gas  temperatures  in 
the  next  increment  are  found  using  equations  (4)  and  (S),  and  the  process 
is  repeated.  A  discussion  of  the  convection  and  radiation  coefficients, 
and  assumption  regarding  them,  is  in  order. 

Convection  heat-transfer  coefficients.  -  The  convection  coefficients 
are  detennined  from  one  of  two  relations  for  Nusselt  number,  depending 
upon  whether  the  flow  in  a  particular  passage  is  laminar  or  t\u:bulent. 

The  relation  for  laminar  flow  is  obtained  from  reference  1: 

^  =  3.  75 
kf 

For  parallel  plates  of  large  aspect  ratio; 

D  w  2s 


Thus, 


7.50^ 


(6) 


For  turbulent  flow  the  equation  for  the  modified  Boelter  equation  given 
in  reference  2  is  used; 


8 


kf 


which  upon  rearrangement  becomes,  for  gases: 


h  - 


0.023 


dO.  2 
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(7) 


In  both  equations  (6)  and  (?)  the  subscript  f  denotes  film  conditions. 
That  is,  the  properties  are  evaluated  at  temperature  tf,  halfway  be¬ 
tween  the  bulk  temperature  of  the  gas  and  the  surface  temperature  of  the 
plate  I 


^  2 

The  validity  of  equation  (?)  has  been  experimentally  verified  for  high 
surface  temperatures  (2000°  to  5000°  R)  and  high  surface  to  gas  tempera¬ 
ture  ratio  (from  2  to  4)  in  reference  3.  Unfortmately,  similar  verifi¬ 
cation  of  the  film  correlation  does  not  exist  for  laminar  flow  (eq.  (6)). 
As  stated  earlier,  the  use  of  equation  (6)  or  (?)  depends  upon  whether 
the  flow  in  a  passage  at  a  particuleo*  increment  is  laminar  or  turbulent. 
It  is  entirely  possible  that  the  flow  in  a  given  passage  cein  change  from 
laminar  to  turbulent,  or  vice  versa,  along  the  length  of  the  passage. 

In  readity  such  transitions  in  flow  are  usually  neither  well  defined  nor 
smooth.  To  avoid  analytical  and  computational  difficulties  associated 
with  determining  the  actual  flow  transition,  a  model  has  been  postulated, 
whereby  the  Nusselt  number  varies  continuously,  as  a  function  of  Reynolds 
number  from  laminar  flow,  through  transition,  to  turbulent  flow.  This 
model  is  illustrated  in  figure  3(a). 

Another  assianptlon  made  regarding  convection  heat  transfer  is  that, 
while  equations  (6)  and  (?)  apply  strictly  to  the  case  where  the  surface 
bounding  the  flow  is  at  a  imiform  temperature,  the  results  may  adso  be 
applied  to  the  present  case  where  the  two  plates  bounding  any  passage  are 
not  generally  at  the  same  temperature. 

The  last  item  of  importance  relating  convection  heat  transfer  con¬ 
cerns  the  entrance  effect  on  heat-transfer  coefficient.  This  effect  re¬ 
sults  in  the  heat-transfer  coefficient  being  increased  over  its  normal 
value  for  approximately  the  first  10  diameters  of  length  along  each  pas¬ 
sage  from  the  inlet  of  the  heater.  This  effect  is  described  for  the  case 
of  tiarbulent  flow  in  reference  4.  For  laminar  flow  the  effect  is  similar 
qualitatively  but  is  somewhat  more  prono\mced.  This  increase  in  heat- 
transfer  coefficient  is  predominantly  a  f\mction  of  the  number  of  diame¬ 
ters,  x/D,  from  the  inlet,  but  it  is  also  a  weak  function  of  Reynolds 
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number.  For  the  present  purpose,  the  dependence  of  the  increase  in  heat- 
transfer  coefficient  on  Reynolds  number  is  neglected,  and  an  empirical 
equation  representing  the  average  results  from  reference  4  has  been  used: 


Equation  (s)  is  assumed  to  be  valid  for  all  Reynolds  numbers,  even  for 
laminar  flow,  and,  although  this  approximation  for  the  case  of  laminar 
flow  is  perhaps  not  very  good,  the  effect  of  this  discrepancy  on  the  over¬ 
all  perfoinance  of  the  heater  is  negligible. 

Radiation  form  factor.  -  The  radiation  form  factor  ^  is  a  func¬ 
tion  of  the  geometry  of  the  enclosure  in  which  the  radiant  energy  is 
transferred  and  of  the  emissivity  and  absorptivity  of  the  radiating  sur¬ 
faces  .  The  first  assumption  that  is  made  in  this  respect  is  that  the 
radiating  surfaces  constitute  "gray  bodies,"  that  is  to  say  that  the 
absorptivity  and  emissivity  of  the  body  are  equal,  which  is  a  good  as¬ 
sumption  for  most  of  the  refractory  metals.  The  second  assumption  is 
that  the  gas  flowing  between  the  plates  is  transparent  to  radiation. 

This  assumption  is  valid  for  gases  such  as  hydrogen  whose  molecules 
are  symmetrical  (p.  82  of  ref.  5). 

For  the  purposes  of  simplification,  it  is  assiuned  that  the  form  fac¬ 
tor  between  two  adjacent  parallel  surfaces  in  an  axial  segment  of  the 
heater  may  be  approximated  by  the  form  factor  for  a  pair  of  infinite  par¬ 
allel  flat  plates.  At  first  glance  this  assimiption  appears  rather  ludi¬ 
crous,  since  the  length  of  each  segment  is  far  from  Infinite j  however, 
on  closer  examination,  this  assumption  is  not  too  bad  if  the  aspect  ratio 
of  each  segment  is  fairly  large  (lO  or  greater),  and  if  the  surface  tem¬ 
peratures  do  not  vary  too  raidically  from  one  axial  segment  to  the  next. 
Aspect  ratio  here  is  defined  as  the  ratio  of  segment  length  to  spacing 
between  the  plates.  The  last  assumption  is  really  an  extension  of  the 
previous  one.  It  is  assimied  that  the  losses  of  radiation  from  the  ends 
of  the  heater  are  negligible.  Reference  6  shows  that  this  is  a  good  as¬ 
sumption,  especially  when  the  temperature  of  the  environment  surrounding 
an  end  of  the  heater  is  nearly  the  same  as  the  surface  temperature  at  the 
end  of  the  heater.  On  the  basis  of  the  infinite  parallel  plate  assump¬ 
tion,  the  form  factor  is  derived  by  McAdams  (p.  63  of  ref.  5)  as 


(9) 
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where  ^n+1  emissivities  of  the  two  adjacent  surfaces  and 

a  is  the  Boltzmann  constant  cr  =  4.  75X10“13  Btu/sec*fM:^-°R^ 

Effect  of  plate  conduction.  -  To  further  simplify  the  solution,  it 
is  assumed  that  the  conduction  of  heat  through  the  plate  material  does 
not  significantly  alter  the  problem.  Specifically,  this  Infers  that  the 
plates  are  sufficiently  thin  so  that  the  resistance  to  heat  flow  across 
the  width  of  each  plate  is  negligible.  On  the  basis  of  this  stipulation, 
the  temperatures  on  both  sides  of  a  plate  may  be  considered  identical. 

In  addition,  if  the  plates  are  thin,  the  resistance  to  heat  flow  along 
the  length  of  each  plate  will  be  very  high,  and  therefore  it  may  be  as¬ 
sumed  that  no  heat  is  conducted  axially  along  the  length  of  the  plate. 

Heat-transfer  summary.  -  The  heat-transfer  analysis  is  now  briefly 
siunmarized.  The  radiation  form  factor  from  equation  (9)  is  substituted 
in  equations  (3).  The  Reynolds  number  in  each  passage  is  computed  at 
each  increment,  and  on  this  basis  the  heat-transfer  coefficients  for  each 
passage  are  computed  from  equation  (6)  or  (?)  in  accordance  with  fig-, 
ure  3.  The  system  of  equations  (3)  is  then  solved  simultaneously  as  dis¬ 
cussed  previously. 


Pressure  Drop 

When  the  heat-transfer  problem  for  the  entire  length  of  the  heater 
has  been  solved,  one  may  proceed  to  determine  the  pressure  drop  through 
each  passage.  Since  the  pressure  drop  in  each  passage  is  dependent  upon 
the  variation  of  gas  temperature  in  that  passage,  it  will  be  necessary 
to  refer  to  the  results  of  the  heat-transfer  calculations  to  obtain  the 
gas  temperatures  upon  which  the  pressure  drops  are  dependent.  As  in  the 
case  of  vhft  heat-transfer  calculations,  the  heater  is  divided  into  100 
axial  segments.  The  pressure  drop  is  then  computed  recursively  over  each 
segment  and  the  total  pressure  drop  in  a  passage  is  then  the  sum  of  the 
pressure  drops  of  all  the  segments  of  that  passage. 

The  pressure  drop  of  gas  flowing  in  a  passage  is  due  to  two  separate 
components.  The  first  component  is  a  friction  component  and  represents 
the  pressure  drop  caused  by  the  shear  stress  between  the  surfaces  bound¬ 
ing  the  passage  and  the  flowing  gas.  It  is  given  by  the  following  rela¬ 
tion; 


=  ir  2f  f  (10) 

where  f  is  the  well-known  Fanning  friction  factor  and  Ax  is  the 
length  of  each  incremental  segment.  The  second  component  contributing 
to  the  pressure  drop  is  the  momentum  component  smd  represents  the  pres¬ 
sure  drop  required  to  cause  an  increase  of  the  fluid  velocity  in  the  pas¬ 
sage.  In  the  heater  under  consideration  the  cross-sectional  area  of  a 
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passage  is  a  constant  over  its  entire  length.  The  reason  the  fluid  ac¬ 
celerates  down  the  passage  causing  a  momentum  pressure  drop  is  that, 
since  it  is  a  gas,  its  density  is  decreased  hy  the  drop  in  pressure  due 
to  friction  and  also  by  the  increase  in  temperature  due  to  heating.  In 
order  to  satisfy  continuity,  for  constant  passage  cross-sectional  area, 
as  the  density  decreases  the  velocity  of  the  gas  must  increase,  resulting 
in  a  pressure  drop  due  to  the  momentum  increase.  The  equation  for  momen¬ 
tum  pressure  drop  Ap^nom  is: 

APmom  =  -  |(Vi+l  -  Vi) 

since 


Combining  equations  (lO)  and  (ll)  to  obtain  the  overall  pressure  drop 
over  the  i'^  increment,  one  obtains 


'n,i 


si 

g 


2f  Ax  ^  /  1  _  1  \ 

(^n,i  *  ^n,  i+l\  \^n,  i+1  •^n,  i  / 

2  / 


(12) 


for  each  passage  n  =  1,  2,  •  •  • ,  Nj  all  increments  i  =»  1,  .  .  . ,  100. 
■Hie  density  p  is  computed  from  the  perfect  gas  relation 


P  = 


Rtc 


(13) 


where  tg  is  the  static  temperature  and  is  related  to  the  total  gas  tem¬ 
perature  obtained  from  the  heat-transfer  calculations  by  the  following 
relation: 


(14) 


8uid  the  Mach  number  M  is  given  by 
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In  solving  equations  (12)  to  (15)  for  the  overall  pressure  drop  in  a  pas¬ 
sage,  it  was  found  advisable,  from  a  computational  standpoint,  to  start 
the  calculations  at  the  heater  exit  with  a  known  value  of  exit  pressure 
(Pn  lOl)  ■*'0  progressively  calculate  the  incremental  pressure  drops 
toward  the  heater  inlet  (Pn,lOO#  Pn  99^  •  '  Pn  l)*  overall  pas¬ 

sage  pressure  drop  in  a  passage  is  ihen  found  from; 

^n  “  Pn, 1  ■  Pn,  101 

t? 

The  reason  for  choosing  this  method  of  computing  the  pressure  drops, 
rather  than  starting  with  a  known  inlet  pressure  and  working  towards  the 
exit,  is  that  in  the  latter  case  the  computed  pressures,  are  continually 
decreasing;  and,  unless  the  pressure  at  the  heater  inlet  is  chosen  with 
extreme  care,  it  is  possible  in  the  process  of  computation  to  generate 
negative  pressures,  clearly  an  unrealistic  and  undesirable  result.  It 
was  fo\md  that  this  condition  was  quite  probable  for  high  values  of  exit 
Mach  numbers,  where  overall  pressure  drop  is  extremely  sensitive  to 
changes  in  inlet  conditions.  This  difficulty  is  clearly  avoided  if  one 
starts  the  computation  at  the  heater  exit  and  computes  the  pressure  drop 
in  the  direction  of  increasing  pressure.  Another  advantage  of  working 
with  a  specified  exit  pressure  is  that  it  allows  one  to  solve  the  problem 
for  a  specified  value  of  exit  Mach  number,  an  important  design  parameter. 

Proceeding  in  the  maimer  described  previously,  equations  (12)  to 
(15)  are  solved  iteratively  for  each  increment  of  passage  length  in  con¬ 
junction  with  equation  (16): 

Pn,l  =  Pn,i+1  -  ^Pn,i  (16) 

The  solution  is  an  iterative  one  since  at  each  increment  only  Pn, i+1 
and  thus  Pn,  i+l  known.  The  density  Pn,  i^  which  appears  on  the 

right  side  of  (12)  and  is  a  function  of  Pn, i,  is  unknown.  Briefly,  a 
guess  is  made  at  Pn,  i^  basis  of  which  Pn,  1  is  calculated  from 

(13);  this  value  of  pn,i  is  then  used  in  (12)  to  compute  Apn^i.  Using 
this  value  of  Apj^^j^  a  new  value  of  Pn^i  is  computed  from  (16),  and 
the  iteration  continues  imtil  the  changes  in  the  computed  vadues  of 
Apn, i  become  arbitrarily  small. 

The  overall  pressure  drops  in  each  of  the  passages  were  computed 
individually  in  the  previous  manner.  The  last  item  to  be  discussed  in 
this  area  is  the  calculation  of  the  friction  factor  f  appearing  in 
equations  (lO)  and  (12).  As  in  the  case  of  heat-transfer  coefficient, 
the  friction  factor  was  also  assumed  to  be  a  function  of  the  film  Reyn¬ 
olds  number.  For  the  case  of  turbulent  flow,  this  ass\miption  is  verified 
by  existing  experimental  evidence  (ref.  2).  Again  the  validity  of  this 
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assumption  for  laminar  flow  is  not  as  clearly  defined.  For  laminar  flow 
the  friction  factor  for  flow  between  infinite  parallel  plates  is 
(eq.  4-85  of  ref.  7); 


f  = 


24 

Re^f 


(17) 


Though  the  parallel  plates  of  this  heater  are  by  no  means  infinite,  the 
aspect  ratio  of  the  cross  section  of  a  heater  passage  is  assumed  large 
enough  that  equation  ( 17 )  is  a  good  approximation . 


For  turbulent  flow  a  popular  equation  for  friction  factor  is  the 
simple  Blasius  relation  (eq.  7-64  of  ref.  7): 

0.046 
^  ReO-2 

Unfortxinately,  this  equation  is  only  approximate  over  the  range  of  Reyn¬ 
olds  nimibers  10^^  <  Re  <  2X10^.  For  Re  <  10^^  the  approximation  of  the 
previous  equation  to  experimental  results  is  particularly  bad. 

A  more  exact  relation  between  friction  factor  and  Reynolds  number 
for  circular  pipes  is  the  so-called  Von  Karman  -  Nikuradse  equation  (eq. 
7-60,  ref.  7): 


=  4.0  logio  (He>^)  -  0*40 
vf 

A  similar  relation  valid  for  flat  plates  may  be  obtained  using  the  same 
methods  (p.  169,  ref.  7)  as  those  used  in  the  derivation  for  the  circu¬ 
lar  pipe  equation  above.  The  Von  Karman  -  Nikuradse  equation  for  flat 
plates  becomes; 

i  =  4. 0  logn o(ReVT)  -  0.  98 
yf 

Obviously,  because  the  friction  factor  appears  on  both  sides  of  the  Von 
Karman  -  Nikuradse  equation,  this  relation  is  not  generally  a  convenient 
way  of  obtaining  friction  factor  as  a  function  of  Reynolds  number.  There¬ 
fore,  it  was  decided  to  develop  an  empirical  equation  that  would  closely 
approximate  the  Von  Karman  -  Nikuradse  equation  for  flat  plates  and  would 
give  friction  factor  explicitly  as  a  function  of  Reynolds  number.  Such 
a  relation  is  given  by  the  following  equation: 

^  ^  _ 0.84 _ 

(in  Ref.  -  1.11)2*22 

Equation  (18)  yields  results  substantially  more  accurate  than  those  given 
by  the  Blasius  relation,  particularly  at  low  turbulent  Reynolds  nimibers. 


14 


As  in  the  case  of  heat  transfer^  to  avoid  computational  difficulties  it 
was  assumed  that  the  variation  in  friction  factor  from  the  laminar  case 
(eq.  (17))  to  turbulent  (eq.  (18))  was  continuous  as  shown  in  figure  3(b). 
As  in  the  case  of  heat-transfer  coefficient,  the  friction  factor  is  com¬ 
puted  on  the  basis  of  the  local  film  Reynolds  number  computed  at  each  seg¬ 
ment  of  each  passage. 


Pressure  Drop  Equalization 

In  obtaining  the  solution  to  the  heat-transfer  and  pressure  drop  prob¬ 
lems,  if  the  mass  flow  is  initially  assumed  to  be  uniformly  distributed 
between  the  passages,  the  resulting  pressure  drops  in  all  the  passages 
will  not  generally  be  equal.  This  is  the  case  even  when  the  dimensions 
of  the  passages  are  identical  and  the  heat  generation  in  each  plate  is 
the  same,  because  the  outer  passages  adjacent  to  the  walls  bounding  the 
heater  element  effectively  receive  heat  from  only  one  side,  since  the 
wall  is  usually  adiabatic  or  at  least  nearly  so.  Since  the  heat  balance 
equations  for  each  plate  are  coupled,  the  fact  that  the  heat  transfer  to 
the  outer  passages  is  not  the  same  as  that  to  the  remaining  passages  af¬ 
fects  the  heat  transfer  and  resulting  temperature  distribution  in  every 
passage  of  the  heater.  This  then  results  in  unequal  pressure  drops  in 
the  passages  despite  the  apparent  symmetry  of  mass  flow,  dimensions,  and 
so  forth  between  the  passages. 

However,  it  is  expected  that,  in  actual  operation,  the  passage  mass 
flows  will  be  distributed  such  that  the  pressure  drops  in  all  the  passages 
will  be  Identical.  The  present  problem  then  is  to  determine  a  distri¬ 
bution  between  the  passages  of  mass  flow  such  that  the  resulting  pres¬ 
sure  drops  will  be  equal.  The  pressure  drop  in  a  given  passage  is  pri¬ 
marily  a  function  of  the  mass  flow  through  that  passage;  because  of  the 
coupling  between  the  passages,  it  is  also  a  weak  function  of  the  mass 
flows  in  other  passages,  but  this  effect  is  second  order  and  is  neglected 
in  order  to  simplify  the  analysis.  Thus, 


Pn  “  P('^n) 


where  is  the  overall  passage  pressure  drop  and  w^  is  the  flow  rate 

in  that  passage,  and  thus. 


where  S^n/^n  influence  coefficients  defined  in  equation  (23). 
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The  APn  represent  the  changes  in  pressure  drop  for  each  passage  re¬ 
quired  to  equalize  all  the  passage  pressure  drops,  and  the  Aw^  are  the 
estimated  changes  in  passage  weight  flows  required  to  yield  the  necessary 
values  of  APn*  At  the  end  of  each  sequence  of  heat  transfer  and  pres¬ 
sure  drop  calculations,  the  AP^  required  to  equalize  the  P^  are  com¬ 
puted.  From  this  one  obtains  the  required  values  of  Aw^  from  equation 
(19),  and  new  values  of  passage  weight  flow  are  computed  by: 

Using  these  new  values  of  passage  weight  flows,  the  heat-transfer  and 
pressure  drop  problems  are  again  solved,  and  this  iterative  procedure 
continues  mtil  the  differences  among  any  of  the  passage  pressure  drops 
become  less  than  some  arbitrary  limit. 

Since  it  is  desired  to  correct  all  the  passage  pressure  drops  to  the 
same  value  F,  AP  in  equation  (19)  is  defined  as 

APn  =  P  -  Pn  (20) 

where  P  is  a  mean  pressure  drop  to  which  it  is  desired  to  correct  the 
Pn-  The  next  consideration  is  by  what  criterion  should  one  compute  P? 

At  first  glance  a  simple  arithmetic  average  of  the  Pn  would  appear  ade¬ 
quate;  however,  this  is  not  convenient  because  of  the  following  addition¬ 
al  consideration.  In  finding  the  solution  to  the  overall  heater  problem 
for  a  particular  case,  the  mass  flow  through  the  entire  heater  is  given 
either  directly  or  indirectly,  and  naturally  this  total  mass  flow  must 
remain  fixed  during  all  phases  of  solution  of  the  problem.  This  imposes 
a  constraint  on  our  pressure  drop  equalization  scheme.  The  perturbation 
of  the  passage  weight  flows  to  equalize  the  passsige  pressure  drops  is 
subject  to  the  constraint  that  the  sum  of  the  passage  weight  flows  re¬ 
mains  constant.  This  can  also  be  expressed  by  saying  that  the  sum  of 
perturbations  in  passage  weight  flows  must  be  equal  to  zero.  Thus, 


=  0 


Substituting  equations  (19)  and  (20)  into  (21)  gives 


(21) 
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(22) 


The  influence  coefficients  used  in  equation  (22)  are  obtained  from  numer¬ 
ical  results.  That  is, 


where  j  designates  the  number  of  the  trial.  Originally,  an  arbitrary 
perturbation  is  made  in  the  values  then  both  the  heat  transfer  and 

pressure  drop  problems  are  recomputed,  and  the  nev  set  o'f  Pj,  values  is 
obtained.  Then  the  first  set  of  influence  coefficients  (Bwn/bPn)]^  is 

computed.  For  the  second  and  subsequent  trials  the  Awj^  are  computed 
from  equation  (22),  the  problem  is  recomputed  using  the  new  mass  flows, 
and  this  process  is  continued  until  the  maximum  difference  between  any 
of  the  pressure  drops  becomes  arbitrarily  small.  The  problem  is  then 
considered  completed.  It  is  important  to  note  that,  in  the  process  of 
iteration  to  equalize  the  pressure  drops,  both  the  heat-transfer  and 
pressure  drop  problems  are  recomputed,  not  just  pressure  drop  alone. 


RESULTS  AND  DISCUSSION 

Results  are  presented  for  the  solution  of  the  heat-transfer  and 
pressure  drop  equations  of  the  heater  using  hydrogen  as  the  gas  with  a 
fixed  heater  inlet  temperature  of  200°  R  and  a  mixed-mean  exit  tempera 
ture  of  5000°  R.  It  was  necessary  to  define  a  mixed-mean  exit  tempera¬ 
ture  since  the  exit  gas  temperatures  in  the  individual  passages  were  not 
generally  equal.  The  hydrogen  pro:perties  used  in  the  computation  were 
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taken  from  reference  8.  Incipient  dissociation  of  the  hydrogen  caused 
an  attendant  sharp  increase  in  specific  heat  and  thermal  conductivity  at 
temperatures  above  4000°  R.  The  effect  of  this  on  the  results  is  to  in¬ 
crease  the  coefficient  of  heat  transfer  near  the  exit  of  the  heater.  At 
temperatures  up  to  5000°  R,  the  effect  of  dissociation  on  the  density  of 
the  dissociated  hydrogen  mixture  is  negligible.  The  results  presented 
here  v/ere  computed  under  the  additional  assumptions  that  the  power  gener¬ 
ation  in  the  plates  varied  axially  according  to  a  2/3-cosine  power  dis¬ 
tribution  (described  in  fig.  4),  that  power  generation  was  equally  di¬ 
vided  among  the  plates,  that  the  spaclngs  between  the  plates  were  identi¬ 
cal,  and  that  the  outer  wall  bounding  the  plate  bundle  was  adiabatic. 

In  addition,  most  of  the  results  are  for  the  case  of  five  plates,  al¬ 
though  results  are  presented  for  nine  plates  for  the  purpose  of  illus¬ 
trating  the  effect  of  changing  the  number  of  plates  on  heater  performsince. 

It  was  felt  that  the  most  significant  parameters  resulting  from  the 
computations  were  maximum  plate  temperature  and  overall  pressure  drop. 
Maximum  plate  temperature  was  selected  as  a  parameter  because,  at  the 
high  operating  temperatures  under  consideration,  the  performance  of  the 
heater  is  limited  by  the  maximuuti  allowable  temperat\ire  for  the  materials 
from  which  the  plates  are  fabricated.  Thus,  the  effect  of  number  of 
plates,  plate  emissivity,  and  L/D  ratio  on  the  maximum  plate  tempera¬ 
ture  was  investigated.  The  results  are  indicated  on  plots  of  maximum 
temperature  against  (L/D)/(GD)®* ^  for  turbulent  flow  and  (L/d)/GD  for 
laminar  flow.  These  two  parameters  are  proportional  to  f(L/D)  for 
tiirbulent  and  laminar  flow,  respectively,  where  f  is  the  txirbulent  or 
laminar  friction  factor,  D  is  the  equivalent  diameter  of  all  the  equally 
spaced  passages,  and  G  is  the  mean  value,  over  the  passages,  of  the 
mass  flow  through  the  heater.  Thus,  one  would  expect  pressure  drop  to 
be  a  function  of  these  parameters,  and,  because  of  the  Reynolds  analogy, 
one  might  also  expect  the  heat-transfer  coefficient  (and  thus  the  maxi¬ 
mum  plate  temperatures)  to  also  be  a  function  of  these  parameters.  Be¬ 
cause  of  the  different  flow  conditions,  the  results  for  laminar  flow  are 
presented  separately.  It  was  felt  that  the  most  probable  region  of  oper¬ 
ation  of  this  heater  would  be  turbulent  flow,  and  therefore  most  of  the 
results  presented  here  are  for  that  case.  Ihe  results  presented  in  this 
report  are  of  course  only  valid  for  the  case  of  hydrogen  as  the  working 
fluid. 

All  the  results  for  maximum  plate  temperature  are  derived  from  plots 
of  the  form  of  figure  5,  which  show  the  variation  of  center  plate  temper¬ 
ature  and  adiabatic  wall  temperature  with  dimensionless  distance  from  the 
inlet.  As  expected,  the  maximum  temperature  for  the  entire  heater  always 
occurred  on  the  center  plate;  since  the  power  generation  in  each  plate 
was  the  same,  the  passages  were  of  equal  cross  section,  and  the  outer 
wall  was  adiabatic.  In  most  cases  the  maximum  plate  temperature  occurred 
in  the  vicinity  of  the  heater  exit.  For  certain  cases  the  meocimim)  point 
was  inward  from  the  exit  because  of  the  combined  effect  of  reduced  heat 
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flux  at  the  exit  given  by  the  2/3-cosine  axial  power  distribution  and  the 
increased  heat-transfer  coefficient  due  to  incipient  dissociation  of  the 
high-temperature  hydrogen.  The  maximum  temperatures  from  plots  of  tem¬ 
perature  against  length  such  as  figure  5  were  obtained  for  a  wide  range 
of  values  of  l/D  and  GD  and  are  plotted  in  the  form  of  figure  6  for 
turbulent  flow  and  in  figure  11  for  laminar  flow. 

In  figure  6  it  is  seen  that  the  parameter  (L/D)/(GD)®* which  is 
proportional  to  f(L/D)  for  turbulent  flow,  does  indeed  represent  the 
heat-transfer  results  for  turbulent  flow  quite  well.  There  are  a  few 
points,  however,  that  warrant  discussion.  It  is  seen  that  the  maximum 
adiabatic  wall  temperature  decreases  as  the  meocimum  plate  temperature 
increases  and  is  a  function  not  only  of  f(L/D),  but  also  of  mass  flow 
G.  The  reason  for  this,  as  will  be  discussed  in  more  detail  later,  is 
that  the  flow  distribution  between  the  passages  is  a  function  of  the 
total  average  mass  flow  G.  As  the  mass  flow  is  increased,  more  flow  is 
shimted  to  the  passsige  adjacent  to  the  wall,  resulting  in  a  drop  of  max¬ 
imum  temperature  there.  Another  interesting  point  is  the  sharp  corner 
in  the  curve  of  maximum  plate  temperature.  This  comer  results  because 
the  absolute  maximum  plate  temperature  shifts  from  a  region  near  the  exit 
to  the  vicinity  of  the  inlet  of  the  heater.  That  the  maximum  plate  tem¬ 
perature  should  ever  occur  in  the  vicinity  of  the  inlet  of  the  heater, 
where  the  cold  gas  is  entering,  does  at  first  glance  seem  surprising. 

This  phenomenon  will  be  the  subject  of  the  next  paragraph. 

The  reason  the  plate  temperature  near  the  inlet  of  the  heater  may 
rise  enough  to  become  the  highest  in  the  entire  heater  may  be  understood 
by  examining  equation  (?).  The  right  side  of  (?)  contains  the  factor 
(tb/tf)®*®.  From  the  definition  of  tf  it  is  seen  that  this  factor  be¬ 
comes  quite  small  when  the  fluid  bulk  temperature  is  low  and  the  surface 
temperature  is  high.  Because  of  this,  the  value  of  h  computed  from 
equation  (?)  is  decreased.  From  equation  (2)  it  is  also  seen  that 
q  =  h(T  -  t^j),  so  that,  as  the  heat-transfer  coefficient  is  decreased, 
the  surface  temperature  must  increase  to  transfer  the  same  heat  flux. 

This  effect  is  also  somewhat  self-perpetuating  since,  as  the  surface 
temperature  is  increased,  (t-b/tf)®’®  decreases  even  more,  resulting  in 
further  reduction  of  heat-transfer  coefficient.  This  effect  is  only  pro¬ 
nounced,  however,  at  high  heat  fluxes  and/or  low  inlet  temperatures,  which 
result  in  low  values  of  t^/tf,  on  the  order  of  0.5  or  less. 

Returning  to  figure  6,  it  is  seen  that  these  high  plate  temperatures 
in  the  inlet  region  occur  at  comparatively  low  values  of  (L^)/(GD)®* 

Since  the  inlet  and  exit  temperatures  for  the  gas  being  heated  sure  fixed, 
low  values  of  l/d  and/or  high  values  of  GD  correspond  to  high  heat 
fluxes  resulting  in  the  phenomenon  described  above.  In  figure  7  one  may 
see  directly  the  effect  of  t^/tf  caused  by  a  change  in  inlet  tempera¬ 
ture  on  the  msucimum  plate  temperature.  The  two  cases  illustrated  are 
identical  except  for  inlet  temperature,  and  the  results  away  from  the 
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inlet  are  nearly  the  sane.  However,  near  the  inlet,  the  plate  tempera¬ 
ture  is  much  higher  for  the  case  of  lower  inlet  temperature.  However,  it 
should  he  mentioned  that  this  effect  at  very  low  values  of  %/tf  has 
never  received  experimental  verification. 

Figure  8  shows  the  effect  of  the  parameter  (L/d)/(GD)®' ^  and  mass 
flow  on  the  flow  distribution  between  the  plates,  which  is  represented  by 
the  ratio  of  the  flow  in  the  passage  adjacent  to  the  adiabatic  wall  to 
that  in  the  innermost  passage.  It  is  seen  that  the  effect  of  mass  flow 
is  only  important  at  conditions  corresponding  to  low  values  of  f(L/D). 
Also,  the  weight  flow  ratio  increases  with  decreasing  f(L/D)  and  increas¬ 
ing  mass  flow.  These  results  are  consistent  with  the  variation  of  adia¬ 
batic  wall  temperature  with  mass  flow  and  f(L/D)  illustrated  in  figure  6. 
The  reasons  for  this  nonuniform  distribution  of  flow  over  the  passages 
are  interesting  and  warrant  some  discussion.  Since  the  outer  passage, 
being  heated  on  only  one  side,  runs  cooler  than  the  inner  passage,  there 
is  less  pressure  drop  there  for  the  same  flow  than  in  the  inner  passages. 
As  a  result,  more  flow  is  shifted  towards  the  outer  passage  from  the  oth¬ 
ers  to  equalize  the  pressure  drop.  As  this  is  done,  the  inner  passages 
run  even  hotter  because  of  their  reduced  flow.  The  situation  is  finally 
equalized  by  the  heat  transfer  from  the  hot  inner  plates  across  the  heat¬ 
er  to  the  cooler  outer  wall  by  means  of  both  convection  and  radiation. 

The  increase  of  the  passage  weight  flow  ratio  with  decreasing 
(lA))/(GD)^* 2  is  believed  to  be  due  to  the  higher  plate  temperatures  as¬ 
sociated  with  low  values  of  (L/D)/(GD)®’ 2.  As  the  plate  temperatures  are 
increased,  the  exit  gas  temperatures  in  the  inner  passages  also  rise,  in- 
c  .'easing  the  pressure  drops  relative  to  those  in  the  outer  passages. 

Thus,  more  flow  is  shifted  from  the  inner  to  the  outer  passage  to  again 
equalize  the  pressure  drops.  The  reason  that  the  passage  flow  ratio  adso 
rises  with  increased  mass  flow  appears  to  be  caused  by  the  fact  that, 
since  the  gas  temperature  at  the  exit  of  the  outer  passage  is  much  less 
than  at  the  exit  of  the  inner  passages,  the  rate  of  increase  of  pressure 
drop  with  increasing  flow  rate  is  less  for  the  outer  passage  than  for  the 
others.  Therefore,  as  additional  flow  is  passed  through  the  heater,  pro¬ 
portionately  more  and  more  of  it  is  passed  through  the  outer  peissage. 


Effect  of  Plate  Emissivity  and  Number  of  Plates 

The  effect  of  changing  plate  emissivity  on  the  maximum  temperature 
is  seen  in  figure  9.  The  results  indicate  that  as  emissivity  is  de¬ 
creased  the  maxlmiim  plate  temperature  is  Increased,  while  the  meoclmum 
wall  temperature  decreases.  Iil  short,  the  temperature  difference  between 
the  center  plate  and  the  wall  is  Increased.  The  reason  for  this  is  that 
decreasing  the  emissivity  increases  the  resistance  to  radiation  heat 
transfer  between  the  plates,  and  thus  the  temperature  distribution  from 
the  center  plate  to  the  waO.!,  across  the  heater,  becomes  more  nonviniform. 
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Increasing  the  number  of  plates  has  an  effect  similar  to  decreasing 
the  plate  emissivity,  in  that  resistance  to  heat  flow  across  the  plates 
is  increased;  this  is  seen  in  figure  10.  The  important  difference  here 
is  that,  while  causing  an  increase  in  the  temperature  difference  between 
the  center  plate  and  the  wall,  increasing  the  number  of  plates  also 
causes  the  maximum  plate  temperature  to  drop.  This  is  in  contrast  to  the 
results  obtained  when  the  emissivity  was  decreased,  which  resulted  in  an 
increase  in  maximum  plate  temperature.  The  reason  for  this  latter  effect 
is  that,  with  a  larger  nxjmber  of  plates  at  the  same  plate  spacing,  a 
larger  percentage  of  the  total  flow  passes  through  the  passages  between 
the  interior  plates,  mitigating  the  temperature  increases  caused  by  the 
passages  along  the  outer  walls  taking  flow  away  from  the  interior  pas¬ 
sages.  As  a  result  the  maximum  plate  temperatures  are  decreased. 


Effect  of  Laminar  Flow 

For  purpose  of  this  discussion,  laminar  flow  operation  is  arbitrar¬ 
ily  defined  to  exist  whenever  the  flow  in  any  of  the  heater  passages  be¬ 
comes  laminar.  Therefore  the  case  of  mixed  flow  conditions  for  the  dif¬ 
ferent  passages  is  included  in  this  definition. 

The  results  for  laminar  flow  are  presented  as  a  function  of  the  pa¬ 
rameter  (l/D)/(GD),  which,  as  for  (l/D)/(GD)®* 2  for  turbulent  flow,  cor¬ 
responds  to  f(L/D).  The  plot  of  maximum  temperature  is  seen  in  figure 
11.  Since  in  the  case  of  laminar  flow  the  heat  fluxes  transferred  in  the 
heater  are  much  lower  than  in  the  turbulent  case,  the  region  of  maximum 
plate  temperature  remains  in  the  vicinity  of  the  heater  exit,  and  no 
abrupt  shifts  of  maximum  temperature  to  the  region  near  the  inlet  of  the 
heater  occur,  as  in  the  turbulent  case.  Another  point  of  interest  is 
that,  except  for  very  low  values  of  GD  (corresponding  to  distinctly  lam¬ 
inar  flow),  the  variation  of  maximum  temperature  with  the  parameter 
(L/D)/GD  is  not  as  well  defined  as  in  the  turbulent  case.  The  primary 
reason  for  this  is  that  at  higher  values  of  GD  the  flow  through  the 
heater  is  no  longer  purely  laminar.  In  figure  12  it  is  seen  that  the 
relatively  large  amount  of  flow  is  shifted  towards  the  outer  passage  in 
the  laminar  cases,  with  the  result  that,  while  the  flow  in  the  inner  pas¬ 
sages  may  be  purely  laminar,  the  flow  in  the  outer  passage  may  be  turbu¬ 
lent.  With  the  heater  operating  in  this  hybrid  flow  regime  the  parameter 
(L/D)/GD  no  longer  can  be  used  satisfactorily  to  represent  the  operating 
characteristics  of  the  heater,  and  the  divergence  of  the  curves  in  fig¬ 
ures  11  and  12  for  low  values  of  (L^)/GD'  is  the  result. 


Overall  Pressure  Drop 

For  the  case  of  turbulent  flow  the  pressure  drop  results  are  illus¬ 
trated  in  figure  13.  Similarly,  the  results  for  laminar  flow  are  plotted 
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in  figure  14.  Again  the  parameters  (L/D)/(GD)®* ^  and  (l/D)/GD^  respec¬ 
tively,  were  used  successfully  in  presenting  the  results.  It  was  found 
that  the  results  were  hest  presented  in  the  manner  shown,  the  pressiire 
drop  being  divided  by  M^p^,  and  the  exit  Mach  number  also  being 

used  as  a  parameter.  The  general  trend  of  results  for  both  laminar  and 
turbulent  flow  is  quite  similar.  In  both  cases  the  pressure  drop  param¬ 
eter  increases  nearly  linearly  with  increasing  f (l/D)  at  constant  exit 
Mach  number  and  the  slope  of  the  pressure  parameter  curves  decreases  with 
increasing  Mach  number.  It  should  be  pointed  out  that,  for  the  cases  pre- 

sented,  MqPq  is  a  parameter  that  is  almost  directly  proportional  to  mo¬ 
mentum  pressure  drop. 

Momentum  pressure  drop  is  given  by  equation  (lla) : 


-  p^)  (lla) 

For  the  cases  presented,  since  Iq  =  5000°  R,  tj^^  =  200°  R,  Po  <  Pin> 
therefore  I/Pq  »  ^/Pin*  Thus  equation  (lla)  may  be  written  approximately: 


■•mom 


-  = 
g  Po 


^ 

g  Po 


(24) 


By  definition. 


Substituting  equation  (25)  into  (24)  one  obtains 

^mom  “  ■  ^PoT 

Otherwise  the  results  for  pressure  drop  are  not  unusual. 


(25) 


SUMMARY  OF  RESULTS 

The  following  results  were  obtained  from  the  analytical  investiga¬ 
tion  of  a  high- temperature  equally  spaced  parallel  flat  plate  heater; 

1.  Heat  transfer  and  pressure  drop  may  be  presented  as  functions  of 
a  parameter  corresponding  to  f(L/D)  for  the  cases  of  both  laminar  and 
turbulent  flow,  where  f  is  the  Panning  friction  factor  and  l/D  is  the 
length-to-equlvalent-diameter  ratio  of  the  passages  in  the  heater. 
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2.  When  the  outer  wall  hounding  the  plate  bundle  is  adiabatic,  the 
outer  passage  receives  more  flow  than  the  interior  passages  in  order  to 
equalize  the  pressure  drops  in  all  passages .  The  proportion  of  the  flow 
passing  through  the  outer  passage  increases  as  f(L/D)  decreases. 

3.  As  the  parameter  proportional  to  f(L/D)  decreases,  the  maxiimmi 
plate  temperature  increases  as  expected.  However,  at  the  same  time  the 
maximum  adiabatic  wall  temperature  decreases  because  of  result  2. 

t 

4.  The  analysis  indicates  that,  at  high  heat  fluxes  (low  values  of  ^ 

the  parameters  corresponding  to  f(L/D))  and/or  low  inlet  gas  temperatures,  c 

it  is  possible  -for  the  maximum  plate  temperature  to  occur  at  a  point  near 

the  inlet  of  the  heater  because  of  the  effect  of  the  temperature  ratio 
t-jj/tf  appearing  in  the  expression  for  the  turbulent-flow  heat-transfer 
coefficient.  However,  this  result  has  never  been  experimentally  verified. 

5.  Increasing  the  number  of  plates  or  decreasing  the  emissivity  of 
the  plates  has  the  effect  of  increasing  the  resistance  to  the  transverse 
flow  of  heat  from  the  inner  plates  to  the  outer  wall,  resulting  in  an 
increase  in  the  difference  between  the  maximum  plate  and  maximum  wall 
temperatures.  As  a  corollary,  the  effect  of  radiation  between  the  plates 
is  to  equalize  the  plate  temperatures. 

6.  Increasing  the  number  of  plates  lowers  the  maximum  plate  tempera¬ 
tures  slightly  because  the  proportion  of  the  flow  passing  through  the 
interior  passages  is  increased. 

7.  The  overall  pressure  drop  through  the  heater  increases  nearly 

linearly  with  increasing  f(L/D)  at  constant  exit  Mach  number  and  pres¬ 
sure,  and  for  a  given  f(L/D)  and  is  roughly  proportional  to  the 

approximate  momentum  pressure  drop  M^Pq- 
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APPENDIX  -  SOLUTION  OF  PLATE  HEAT  BALANCE 
EQUATIONS  FOR  PLATE  TEMPERATURE 


The  set  of  equations  for  the  plate  temperatures  to  he  solved  is: 

'Ign  ~  ^  ^  “  ^n+l)  “  ^-l('^n  "  '*'n-l^  " 

n  =  1,  2,  .  .  .  ,  N  +  1  (3) 

This  system  constitutes  a  set  of  N  +  1  polynomial  equations  in  the 
N  +  1  unknowns  T^^,  T2,  •  •  • ,  Tjj+i*  The  solution  to  this  set  of  equa¬ 
tions  will  he  obtained  using  the  Newton-Raphson  method,  described  in  ref¬ 
erence  9.  An  initial  set  of  trial  values  of  T^  is  substituted  into  the 
set  (3).  Since  these  trial  values  do  not  represent  the  solution  exactly, 
substitution  into  the  left  side  of  (3)  will  result  in  the  right  sides  be¬ 
coming  nonzero.  Designate  the  nonzero  residuals  on  the  right  sides  by 
Rjj.  Then  upon  slight  rearrangement,  (3)  is  rewritten  to  give: 

<lg^  ^n-l'^n-1  *  ^*^-1  ^ n)"^^  ‘  (^-1  ^)'^n  •^^n+l 

+  hj^_i^n-l  ^^n  ~  n  =  l,  2,  •  .  .,N  +  1  (Al) 

The  Newton-Raphson  method  involves  the  solution  of  the  following  set  of 
N  +  1  linear  equations  for  ATjj,  the  changes  to  be  applied  to  the  trial 
values  of  T^  to  give  the  new  trial  values  of  T^.  This  set  of  equa¬ 
tions  is  described  by: 

^n-1  +  ^  +  ^2-  ^n+1  =  -  ^n 

^n-1  '^n+l 

n  =  1,  2,  .  .  N  +  1  (A2) 

The  other  partial  derivatives  SRn/STn-2,  SRn/STn+2^  ®o  forth  do  not 
appear  because  of  the  form  of  equation  (Al) .  Partial  differentiation  of 
(A1)  and  substitution  into  (A2)  give: 

(4^n-lT^-l)^n-l  "  j^C^n-l  + 

-  (Tn  -  tn-l) 

The  terms  of  the  form  dh^/^Tn  ! 
film  temperature  tf^n  =  (tb,n+  Tn)/2,  making  hn  functionally  dependent 
upon  Tn- 


^n)T§  -  (hA-1  +  hn)  '  (^n  -  tn)  ^ 

Win 

ATn+i  =  -  Rn 

ire  due  to  the  dependence  of  on  the 
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As  in  the  main  text,  all  terms  having  subscripts  of  0  or  N  +  2 
in  the  first  or  last  equations,  respectively,  are  identically  zero.  For 
simplicity  the  following  terms  are  now  defined: 

^n  =  4  -^-iT^-l 

hn  =  -  .^n-1  +  ^ n)Tn  “  (^^n-1  +  ^n)  "  ('^n  "  "tn)  ^5^  "  (^n  “  '^n-l) 

'^n  =  n^n+1 

Then  the  system  of  equations  (A3)  is  more  simply  described  by: 

b^  +  cj^  ZYTg  =  -Ri 


®n  ^n-1  +  ^n  ^n  “^n  ^n+1  “  "^^n  n  -  2,  3,  .  .  . ,  N 


=N+1  ^  +  ^N+1  ^N+1  =  ■%+! 


(A4) 


The  solution  of  this  system  of  linear  equations  for  ATn  presents  no 
problem  when  performed  by  matrix  methods  on  a  digital  computer.  However, 
the  use  of  general  matrix  inversion  routines  for  the  solution  K  linear 
simultaneous  equations  on  a  digital  computer  in  this  case,  though  entire¬ 
ly  possible,  is  wasteful  both  of  computer  time  and  storage.  Because  of 
the  particular  arrangement  of  elements  about  the  main  diagonal  of  the  co¬ 
efficient  matrix  for  (A3)  and  the  resulting  large  number  of  zero  elements 
in  that  matrix,  a  direct  elimination  method  may  be  used  to  advantage  over 
the  more  general  technique  in  obtaining  the  solution  to  this  particular 
set  of  equations.  The  direct  elimination  method  may  also  be  programmed 
for  a  digital  computer,  and  it  is  much  simpler,  faster,  and  shorter  than 
a  matrix  inversion  routine.  A  brief  description  of  this  method  has  re¬ 
cently  been  given  in  reference  10. 


When  the  obtained  from  the  solution  of  (A4),  the  new  Tn 

for  the  next  trial  are  computed  as  follows; 


"^n,  Z+1  =  ^n,  Z  ^n,  Z 
n  =  1,  2,  .  .  . ,  N  +  1 
Z  =  Trial  nimiber 


With  the  ne.-r  values  of  T^,  the  calculations  are  repeated,  beginning  with 
equation  (Al),  and  this  process  is  continued  until  all  the  become 

arbitrarily  small,  the  resulting  values  of  Tn  being  a  good  approxima¬ 
tion  to  the  exact  solution  of  (3). 


4 
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Figure  1.  -  Parallel  plate  heat  geometry. 
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Figure  5.  -  Typical  variation  of  plate  and  wall  temperature 
along  length  of  heater.  Axial  power  distribution,  2/3 
cosine;  exit  static  pressure,  40,000  pounds  per  square 
foot;  plate  length,  4.0  feet;  exit  Mach  number,  0.4;  exit 
total  temperature,  5000®  R;  inlet  total  temperatxxre, 

200°  R;  plate  spacing,  0.004  foot;  niamber  of  plates,  five; 
plate  emisslvity,  0.35. 
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Figure  7.  -  Effect  of  high  Tg/t-jj  at  inlet  (low  inlet  tempera¬ 
ture)  on  axial  temperature  distribution  along  plates.  Axial 
power  distribution,  2/3  cosine;  exit  static  pressure,  40,000 
pounds  per  square  foot;  exit  Mach  number,  0. 4;  total  length 
of  plates,  3.0  feet;  plate  spacing,  0.006  foot;  exit  total 
temperature,  5000°  R;  plate  emlssivlty,  0. 35;  number  of  plates 
five. 
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Figure  8.  -  Ratio  of  flow  In  outer  passage  to  that  in  innermost  passage  as  a  function  of  (L/d)/(GD)0* 2  for  turbulent 
flow.  Exit  total  temperature,  5000°  Rj  inlet  total  temperature,  300°  R;  number  of  plates,  five;  plate  emissivlty, 
0.35.  W3  is  flow  in  passage  adjacent  to  adiabatic  wall;  w^  is  flow  in  innermost  passage. 
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Figure  13.  -  Variation  of  pressure  drop  parameter  with  Mach 
number  and  turbulent  flow  correlating  parameter  for  turbulent 
flow.  Inlet  total  temperature,  200°  Rj  exit  total  tempera¬ 
ture,  5000°  R|  nianber  of  plates,  five. 
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Figure  14.  -  Variation  of  pressure  drop  parameter  with  Mach  nianber  and  laminar 
flow  correlating  parameter  for  laminar  flow.  Inlet  total  temperature,  200®  R; 
exit  total  temperature,  5000°  R>  number  of  plates,  five. 
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